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Abstract 

Spin squeezing via atom-field interactions is considered within the context of the Tavis-Cummings 
model. An ensemble of N two-level atoms interacts with a quantized cavity field. For all the atoms 
initially in their ground states, it is shown that spin squeezing of both the atoms and the field 
can be achieved provided the initial state of the cavity field has coherence between number states 
differing by 2. Most of the discussion is restricted to the case of a cavity field initially in a coherent 
state, but initial squeezed states for the field are also discussed. Optimal conditions for obtaining 
squeezing are obtained. An analytic solution is found that is valid in the limit that the number 
of atoms is much greater than unity and is also much larger than the average number of photons, 
a^, inititally in the coherent state of the cavity field. In this limit, the degree of spin squeezing 
increases with increasing a, even though the field more closely resembles a classical field for which 
no spin squeezing could be achieved. 



I. INTRODUCTION 



Spin squeezed states offer an interesting possibility for reducing quantum noise in preci- 
sion measurements [1-3] . Spin squeezing is described in terms of spin operators that are asso- 
ciated with quantum mechanical operators of two-level atoms (TLA) (we refer to atoms and 
spins interchangeably). In an appropriate interaction representation, combinations of atomic 
raising and lowering operators for atom j are associated with the x and y spin components 
{Si and S-j^), while the population difference operator for the two states is associated with the 
z spin component (Si) . One then defines collective operators — S^^ that obey the usual 
spin commutator relations. If one measures an average spin | (S) | = ^ i^x)^ + (Sy)^ + (Sz)'^ 
then the system is said to be spin-squeezed if 



where AS± is the uncertainty in a spin component perpendicular to (S), S — N/2, and N 
is the number of atoms [1, 2]. Spin squeezing is impossible for a single atom and requires 
the entanglement of the spins of two or more atoms. There are many ways to theoretically 
construct a Hamiltonian that can give rise to the necessary entanglement among N two- 
level atoms. Since a linear Hamiltonian merely rotates the spin components leaving the 
uncertainties unchanged, it is generally necessary to use Hamiltonians that are quadratic 
in the spin operators to generate squeezing. On the other hand, it is possible to generate 
squeezing using a Hamiltonian linear in the spin operators provided the spin system is 
coupled to another quantum system, such as a harmonic oscillator. It is then not surprising 
to find that a squeezed state of the oscillator can be transferred to some degree to the atoms. 
What may be a little more surprising is that an oscillator prepared in a coherent state and 
coupled to the spins can result in spin squeezing. In this paper, we study the dynamics 
of the creation of squeezing in an ensemble of spins via coupling to a cavity field in the 
Tavis-Cummings model [4]. An ensemble of N atoms is coupled in a spatially independent 
manner to the N atoms with no losses for the field and with the neglect of any spontaneous 
emission for the atoms. We are concerned mainly with the type of spin squeezing that can 
be generated by coupling to a radiation field that is initially a coherent state, but also will 
consider an initial state of the field that is a squeezed state. The evolution of the radiation 
field will also be determined. There have been a number of studies of atom-field dynamics 
in the Tavis-Cummings model in which the squeezing of the cavity field was calculated in 
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various limits [5]. Some numerical solutions to the problem of spin squeezing in the Tavis- 
Cummings model are given in Ref. [1]. 

The initial condition for the atoms is taken as one in which all the atoms are in their 
lower energy state, corresponding to a coherent spin state. For a very large number of atoms 
(A^ 3> 1 and N much greater than the average number of photons in the coherent state of 
the field) , the relevant energy levels of the spin system approach those of a simple harmonic 
oscillator with corrections that vanish as A?" ~ oo. Thus it would seem that spin squeezing 
can never be achieved if the initial state of the cavity field is a coherent state, since one is 
dealing with a linear interaction between two harmonic oscillators each of them initially in 
a coherent state. Nevertheless, we show that for any finite N, spin squeezing occurs and the 
degree of spin squeezing actually increases with increasing field strength. 

To follow the atom- field dynamics, we consider first a system having N — 2. It is not 
difficult to obtain analytic solutions in this case, enabling us to track the dependence of ^± 
on field strength and N. In addition, we determine if the squeezed vacuum state results in 
optimal transfer of squeezing from the fields to the atoms. After discussing the two atom 
case, we generalize the results to N atoms. 

The paper is organized as follows. In Sec. II we present the mathematical framework 
and obtain results that show that no squeezing can be achieved when the field is either 
classical, or quantized in a number state. In Sec. Ill, we consider the N — 2 case and obtain 
analytical results for both coherent and squeezed cavity fields, in the limit that the average 
number of photons in the field is much less than unity. Numerical solutions for larger field 
strength are presented. In Sec. IV, the results are generalized to A^ atoms. In both sections 
III and IV, the time evolution and squeezing of the field is also calculated for the case that 
the field is initially in a coherent state. In Sec. V, a formal derivation of the large A^ limit is 
given using the Holstein-Primakoff transformation [6], valid for an arbitrary strength of the 
coherent cavity field. The Holstein-Primakoff transformation was used previously by Persico 
and Vetri [7] to analyze the atom-field dynamics in the hmit of large N. The approach we 
follow differs somewhat from theirs and our results seem to have a wider range of validity 
than that stated by Persico and Vetri. The results are summarized in Sec. VI. 
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II. GENERAL CONSIDERATIONS 



In dipole and rotating-wave approximations, the Hamiltonian for an ensemble of TLA 
(lower state |1), upper state |2) , transition frequency cu) interacting with a resonant cavity 
field, E{t) = Ea e"^'^* + E*a ^e''^\ is of the form 

H = hwSz + huia^a + hg{S+a + S^a'^), 
where = E7=i [(|2) (2|), - (|1) /2, = E7=i (|2) (1|), e-*, = 



J2j=i (2|)^- = {S+ + S_) /2, Sy = (S+ - /2i, a and at are annihilation and 

creation operators for the field, and is a coupling constant. The spin operators have been 
defined in a reference frame rotating at the field frequency. Constants of the motion are 
S"^ = 3"^ + Sy + si and {Sz + a^a). If, initially, all spins are in their lower energy state, then 
^2 _ 7\/"2^^ order to calculate from Eq. (1), one must first find (S) and define two in- 
dependent directions orthogonal to (S) , and 5'_l2. It then follows that (-S'j.i) = ('S'_l2) = 
and 

where i = 1, 2 and 5'^-'^ is a spin operator for atom j. 

A necessary condition to have < 1 is that the different spins are entangled. To see 
this, take a system in which (S) is aligned along the z axis, with the x axis is chosen such 
that is the minimum value of S^. Using the facts that (S) = S^i {Sx) = {Sy) — 0, 
ASxASy > \{S,)\/2, one finds 



NASJl (S,) I > VN/ASy 



For correlated states, the sum can be positive and one cannot rule out the possibility that 
< 1. On the other hand, for uncorrelated states, using the fact that (Sy)'^ — 0, it follows 

that 1 + J2j j'^j (^Sy'^^y '*) ^ ~ Ej (^^^y^^ ■ ^ consequence, > ^ and there is no spin 

squeezing for uncorrelated states. 

We note two general conclusions that are valid for arbitrary N. First, if wc were to replace 

the cavity field by a classical field, the Hamiltonian would be transformed into 



where g' is a constant. Since the Hamiltonian is now a sum of Hamiltonians for the individual 
atoms, the wave function is a direct product of the wave functions of the individual atoms. As 
a consequence, there is no entanglement and no spin squeezing for a classical field. Second, 
if the initial state of the field is a Fock state, although there is entanglement between the 
atoms and the field, there is no spin squeezing. There is no spin squeezing unless the initial 
state of the field has coherence between at least two states differing in n by 2. For a Fock 
state, there is no such coherence and ^± > 1. 

It is convenient to carry out the calculations in an interaction representation with the 
wave function expressed as 

N/-2 oo 
m.=-N/2 n=0 

where m labels the value of Sz and n labels the number of photons in the cavity field. In 
this representation, the Hamiltonian governing the time evolution of the c„fc(i) is given by 

H = hg{S+a + S-a+). (3) 



III. N=2 



We first set N — 2, S — 1. If the spins are all in their lower energy state Bit t — 0, the 
initial wave function is 

oo 

wo)) = E^'^i-i'^)' (4) 



k=0 



where the Ck are the initial state amplitudes for the field. Solving the time-dependent 
Schrodinger equation with initial condition (4), one finds 



= k - 1 + A;cos(V4A;-2c/i)] 



k+l 



co,k{t) = -iJ ^j—-^sm{VAkT2gt)ck+i 



Cl,k{t) = 



-1 + cos{V 4:k + 6gt) 



Cfe+2- 



(5a) 
(5b) 
(5c) 



2k + 3 

These state amplitudes can be used to calculate all expectation values of the spin operators. 
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A. Coherent State 



If the initial state of the cavity field is a coherent state, then 

Cfe = a^e-l"l'/VVM, (6) 

and the average number, no, of photons in the field is given by no = |a| . For simplicity, we 
take a and g to be real. 

2 

1. Solution for \a\ « 1 

Keeping terms to order a^, one finds from Eqs. (5) and (6) that the only state amplitudes 
of importance are 



c_i,oW = (1 - «V2) 
c-i,i(i) = a cos {V2gt) 



c-i,2{t)^^[l + 2cos{V6gt) 
Co,o(^) = —i<^sm{^/2gt) 



10^ 

Co,i{t) = --^sin(V6^i) 



ci,o(^) = - Y 1 - cos{\/^gt) 



The spin components' expectation values are: 



{S^) = {Sy) = y/2as\n{V2gt) 



\-Q? sm^{V2gt) 



(7a) 
(7b) 

(7c) 

(7d) 

(7e) 

(7f) 



(8a) 
(8b) 



The motion of the average value for the spin vector operator is in the yz plane, with the 
length of the vector always equal to unity, to order a^. Since (Sx) = 0, the plane in which we 
look for spin squeezing is the one defined by the x axis and an axis orthogonal to both x and 
the instantaneous direction of the spin. Making the appropriate rotation in the yz plane to 
define a y' axes perpendicular to (S) and x, and afterwards choosing an arbitrary direction 
defined by an angle in this plane, one finds that ^ min{ ^j., ^y/}, which implies that 
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FIG. 1: Spin squeezing as a function of gt for a = 0.4 and N = 2. 



the best squeezing is to be found in either the x or y' directions. The analytical expressions 
for are: 



= ^ 1 + { ^ sin2(v^^t) - ^ sin^ [V6gt/2) | (9a) 

e.' = ^ 1 + { - ^ sin^(v^^t) + ^ sin^ [V6gt/2) } (9b) 

The lowest possible value for the squeezing occurs in the x direction and is equal to 

at a time when sin {V2gt) = and cos(V6^i) = -1. The squeezing as a function of gt for 
q; = 0.4 is plotted in Fig. Ill A 1. 



2. Numerical results for all values of a 



General expressions for the spin expectation values and variances can be obtained and 
used for numerical simulations for any values of a. With a real, the expectation value of 
the X component of the spin vanishes and, with the notation cq^ = 



7 



oo 



(Sy) = V2^Co,n(cL 



,n 



C-l,n) 



n=0 



oo 



C-l,n| ) 



n=0 



The variances are: 




(11a) 




} - {Syf 



(lib) 



The variance in the x component of the spin cannot be less than 1/2 unless ci,„c_i,n < 0. 
Since ci,„c_i,n is proportional to Ck+2Ck, where the CfcS are initial state amplitudes for the 
cavity field, spin squeezing can be induced by a field only if the field has at least one 
nonvanishing off-diagonal density matrix element p^k' for which \k — k'\ = 2. 

The values for the spin averages and uncertainties are calculated in terms of a and gt. 
For « 1 the numerical and analytical results agree. For larger values of a, no analytical 
solution is available. The numerical results indicate that the optimal squeezing is obtained 
in the x direction. As a is increased, the spin squeezing increases and then decreases for 
a > 0.9, as shown in Fig. 11. 

With increasing a, the optimal squeezing occurs at increasingly large values of gt. For 
example, with a — 1.6, there is effectively no spin squeezing for gt < 333 and the optimal 
spin squeezing occurs for gt — 2439. The squeezing data in this and subsequent graphs is the 
optimal squeezing that is obtained for gt less than some arbitrary cutoff that we have chosen. 
In the limit of a » 1, the field closely resembles a classical field and {^x)^ia^PP^oach.es 
unity. Formally, this result could be derived by using a transformation proposed by MoUow 
[8] in which the transformed Hamiltonian is that of a classical field having amplitude a 
plus a fluctuating fleld. Any spin squeezing that is produced depends on the ratio of the 
fluctuations to the average fleld strength and must decrease with increasing a, provided the 
average number of photons in the fleld is much larger than N. 
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FIG. 2: Optimal spin squeezing {$,x)min as a function of a for A?^ = 2. The time range out to 
gt = 5000 was explored in obtaining the minimal squeezing. In this and other plots, the point 
represent actual values for which the squeezing was calculated. A line is drawn through these 
points. 




Heisenberg limit 



± 



B. Squeezing in the radiation field 

Although the field is initially in a coherent state, it is squeezed as a result of its interaction 
with the atoms [5]. In terms of quadrature operators P and Q defined as: 

with [Q,P] = i, squeezing of the field occurs if the variance of one of these two operators is 
smaller than the value it would have for the vacuum field. Initially the field is in a coherent 
state of real amplitude a with (^Q^ = \pla and = 0, and variances (AQ)^ = (AP)^ = | 
satisfying the minimum uncertainty condition 



1 



4 



2 1 



(Ag)^(AP)^ = - [g,p] (21) 



4 
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Using the wave function (2), one finds 

1 oo 



m=-l k=0 

1 1 

2 + 2 

1 1 



m=-l A;=0 
1 oo 

' =0 

'5') = i + ^ [^/k{k-l)cl,^k_^{t)Cm,k{t) + ^/{k + l){k + l)4,fc+2(t)c™,fc(t) + 2ks/k + 1 |c^,fc(t)|' 



P')^ 

2 2 



To order a^, for the field initially in a coherent state, one finds squeezing parameters 



= V2AQ ~ 1 - |cos^(V2^t) - - p + 2 cos{V6gt) 
Cp = V2AP ~ 1 + < cos^ {V2gt) - - 1 + 2cos(y65t) 



With this definition, squeezing occurs for < 1 or < 1. To second order in a the 

state of the field evolves in time as a minimum uncertainty state but with squeezing transfer 
between the two quadratures. The minimum value for the squeezing parameters that can 
be obtained is: 



1 - -a' 



(Cp)min — 1 



a 



(12a) 
(12b) 



A continuous transfer of squeezing between the Q quadrature and the x component of the 
spin, and also between the P quadrature and the y component of the spin is taking place. 
The maximum field squeezing as a function of a is shown in Fig. 12. 



C. Squeezed initial cavity field 

^From Eq. (11) one can see that initial state coherence between photon field states 
differing by 2 is needed for squeezing. The squeezed vacuum is a superposition of even Fock 
states; therefore, it is a good choice for inducing the necessary coherences in the atomic 
system. Analytical results are available for a small squeezing parameter of the field, and 
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FIG. 3: Optimal field squeezing {^q)^^^ as a function of a for AT = 2. The time range out to 
gt = 5000 was explored in obtaining the minimal squeezing. 



numerical results can be obtained for larger values. For a squeezed vacuum the c^s are given 

by 

(A; - l)!!(-l)'=/Hanh'=/V 



Co = 



V cosh r ' ^ Vk\ coshr 

where r is the squeezing parameter, assumed real. For any field containing only even expan- 
sion coefficients, (Sx) = (Sx) = 0. For r -C 1, one obtains for the spin squeezing 



for k even; Ck — for k odd. 



^ ASx 4 2 



gt) 
gt) 



To the first order in r, the resulting state is a minimum uncertainty state, and the minimum 
squeezing that can be achieved is the same for both components. Squeezing as a function 
of r is shown in Fig. Ill A 2. 

With increasing r, decreases to minimum value of 0.78 for r ~ 0.7, and then increases 
with increasing r. This result is consistent with the general conclusion that optimal squeezing 
is obtained when the average number of photons in the field is much less than N. 
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Heisenberg limit 
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FIG. 4: Optimal spin squeezing (C?/)min ^ ^ function of the squeezing parameter r for an initially 
squeezed cavity field for N = 2. The time range out to gt = 5000 was explored in obtaining the 
minimal squeezing. 

One might think that the squeezed vacuum produces optimal squeezing, but field states 
that more closely approach the Heisenberg limit = 1/V2 can be constructed. One such 
state is 

1^(0)) = -0.79 |0) - 0.594 |2) + 0.15 |4) + 0.021 16) 

for which a minimum value = 0.724 is achieved. We have not been able to formulate a 
general proof as to the minimum squeezing one can obtain for an arbitrary initial state of 
the field. 

IV. N ATOMS 

As the number of atoms, N, increases, the spin squeezing that can be achieved depends 
critically on the initial state of the cavity field. If the field is in a coherent state, one 
might expect that the squeezing goes to zero as N goes to infinity since the atomic spin 
Hamiltonian approaches that of a simple harmonic oscillator in this limit. A formal proof 
of this result is given below. On the other hand, for finite A^, there are times for which spin 
squeezing occurs, and the squeezing decreases with increasing field strength, provided N is 
much larger than the average number of photons in the field. If the initial state of the field 
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is a squeezed state such as the squeezed vacuum, the field squeezing can be transferred to 
the atoms. In this manner, one can generate a high degree of spin squeezing <^ 1, but 
still considerably less than that predicted by the Heisenberg limit = 1/ y/N . 

For arbitrary A^, the cavity field can, in principle, couple {N + 1) collective states corre- 
sponding to the angular momentum manifold S = N/2. In practice, the number of states 
coupled is on the order of the average number of photons in the initial field. The equations 
of motion for the state amphtudes, obtained from Eqs. (2) and (3) are 




) ~ + ^) + ^) Cm-l,n+l + ~ + m + 1^ n Cm+l,n-l | , 



(13) 

with initial condition Cm,niO) = CnSm-N/2- This equation represents a set of coupled equa- 
tions, starting from m = —N/2 and reaching some maximum value to —N/2 plus n^ax, 
where nmax is the smallest n where the initial field state amplitude c„ is negligibly small. 

A. Coherent cavity field 

1. Analytical solution for \a\^ « 1 

For <S 1, the lowest order non-vanishing amplitudes obtained from Eqs. (13) and (6) 

are 

c_5,o(i) = (1 - «V2) (14a) 
c-s,i{t) ^acos{^/Ngt) (14b) 
c_5+i,o(t) = -iasin(\/iV5t) (14c) 

^-^'^(i) = ^^[N - 1 + Ncos{VW^gt) ] (14d) 

c^s+i,i{t) = -i-^=^= sin(V4Ar - 2gt) (14e) 



y/AN-2 

a^j2N(N-l), , , , 

c_5+2,o(i) = . ^ 1 - cos(V4]V^^i) 



\N -2 

In the large N limit, the average spin components calculated using these amplitudes are 



(14f) 



{S,) = 0; {Sy) = ^fNa%^^{^fNgt)■ (S,) ^-- + a^ sin^VNgt) , (15) 

13 



such that |(S)| = ,5 = N/2 to order a^. 

The squeezing parameter, calculated using Eqs. (11), is given by: 



= + sin^( v^^t) - sin^ (V(2Ar-l)/2,t) } (16) 

In the limit of large N this reduces to 

~ 1 + sin 

As N approaches infinity, the squeezing vanishes; however, for any finite A^, there is a time 
of order 1'k\IN j g where spin squeezing with ~ 1 — occurs. Note that, for small 
gt <S N~^l'^^ ^j. from (16) varies as [l — {gtf{N-l)/&] while from (17) varies as 
\l -\- {gtf /2\^ which have different functional forms; however, the difference between 
these two results varies as /N <^ <^ 1. 

2. Numerical results for all values of a 

Since the average number of photons in a coherent state is a^, one needs to solve Eq. (13) 
up to terms with n ^ cc^. As a grows the numerical solution becomes somewhat unwieldy. 
In Fig. IV A 2, the optimal squeezing is plotted as a function of A^ for a = 0.5. 

The squeezing diminishes with increasing A^, eventually reaching an asymptotic value of 
0.86. This result represents the general trend that the squeezing saturates for A^ 3> a^. Spin 
squeezing as a function of a for fixed A^ = 20 is shown in Fig. IV A 2 

for 0.6 < a < 3.5. The values of in Fig. 6 are do not necessary represent the optimal 
spin squeezing; rather they give first minimum of the envelope of a graph of versus gt. 
It is possible that better spin squeezing occurs at higher values of gt than those considered 
(e.g., for a = 0.6, the first envelope minimum at gt = 9.03 gives = 0.906, while the 
second envelope minimum at gt = 28.1 gives = 0.817); the computation time that would 
be needed to determine (^x)min ^^i all values of gt grows rapidly with increasing a. Spin 
squeezing improves with increasing a up to a ~ 2.7 ~ O(-\/20) and then decreases with 
increasing a, following the general trend noted above.. Spin squeezing for larger values of a 
and N ^ are better treated by the method given in Sec. V. 



2VN- 



gt 



sm 



gt 



(17) 
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FIG. 5: Optimal spin squeezing function of N for a = 0.5. 
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FIG. 6: Optimal spin squeezing (Cx)min ^ function of a for = 20 and < gt < 10. Since only 
a restricted range of gt was considered, the values plotted may not represent the global optimal 
squeezing, but still reflect the qualitative variation of (Cx)min with a. 
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3. Squeezing in the field 



For a <C 1, one finds squeezing parameters 



N - 1 + A^cos(\/4iV-2^t) 
N - 1 + Ncos{-\/AN -2gt) \ , 
implying that 

2N ^ , , 

(eQ)min = 1 - (18a) 
(ep)min = 1 - (18b) 

The best squeezing is obtained ior N — 2. With increasing a, the field squeezing mirrors 
the spin squeezing. 

B. Squeezed initial cavity field 

The spin squeezing one can achieve increases dramatically if the initial state of the cavity 
field is a squeezed state. For a squeezed vacuum with squeezing parameter r, the initial 
squeezing in one quadrature component of the field is — e"'^. In the limit that N ^ 
sinh^r + \/2 sinh r cosh r = (average plus standard deviation of the number of photons in 
the original cavity field), one can show [1] that this squeezing can be transferred totally 
to the spins = <^ '- For large r, this represents substantial squeezing, but since ^ 
sinh^r + sinh r cosh r, it follows that — » ^^0- Thus, one is still far from the 
Heisenberg hmit. It may be possible to construct an original cavity field state that leads 
more closely to the Heisenberg limit — ^/ VN, but we have not explored this possibihty 
in the large N limit. 

V. ASYMPTOTIC SOLUTION FOR LARGE N 

For an ensemble having a number of atoms much larger than unity and much larger than 
the average number of photons in the field, the interaction between the atoms and the cavity 



iQ = V2AQ ~ 1 



2 2/ 

a cos 



Ngt) 



2N -^ 



ip = V2AP ~ 1 + a'^cos'^{VNgt) - — — 
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field can be seen as an interaction between a harmonic oscillator (the field) and an imperfect 
oscillator (the atoms). To attempt to map this problem into one of interacting harmonic 
oscillators, which will be valid as the number of atoms N approaches infinity, one defines 
boson operators for the atoms via 

S, - -N/2 + b^b, (19a) 



(19b) 



S+ = e-'^'N'/^b^{l - b^b/Ny/^ ~ e-^' (^VNb^ - ^^b^b^b^ , 
S- = e^'N'/^l - b^b/Ny/% ~ e^' (^^/Nb - ^^b^bb^ . (19c) 

The boson occupation states (Fock states) \m) = -^^^lO) correspond to the different projec- 
tions onto the collective angular momentum states and, in effect, represent excitations above 
the lowest state having Sz — —N/2. The transformation to the b bosons (Holstein-Primakoff 
transformation [6]) is exact. The approximations in (19a) and (19b) are valid provided the 
relative variations of the spin projection are small: 

{b^b)/N < 1; (20) 

in other words, the average spin remains aligned very close to the z axis. The key point in 
this calculation is that all changes in the eigenkets of order 1/y/N are neglected. Changes 
in the eigenenergies of order 1/\/N lead to significant changes in the phases of the time- 
dependent wave function for any finite A^. Such changes in the phase can result in spin 
squeezing. 

The total Hamiltonian ( in an interaction representation) is written as H — Hq -\- H' , 
with 

= h^fNg {b^a + a^b) , (21a) 
H' = {b^b^ba + a^b^bb) . (21b) 



We now diagonalize Hq and treat H' as a perturbation. The Hamiltonian Hq can be written 
as 
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with 



rt = ^L±l. ^ 

at - 6t _ 

while perturbation if' has the form 

H' = [^t^t^^ _ rtr^rr + r^r^-f-f + -f^-f^rr + 2 {r^-f^i^y-f - rr) + /i.e.}] . 

4v A'^ 

Only the first two terms in this expression, conserving the total number of excitations, 
contribute in first order. 
The eigenkets of Hq are 

\n)^\m). = ^|0)^|0), (22) 



with energies 

e*'^^ (n, m) — h (ou+n + uj-m) . 
The first order correction to the energies of these states is 

, ^ , 2 2\ (1)/ \ (1)/ \ 

'[n^ m) = — -^{n — n — m + m ) = el (nj + el (m). 



and we define 

e±(n) — tujj±n + e^^{n). 

To this order the states (22) are unmodified. 

In order to neglect higher order correction to the energies, it is necessary that the 
phase produced by such corrections must be much less than unity. This translates into 
the condition — -^^j2 ^ 1; which can always be satisfied for sufficiently large 

N, but would be violated ior N — 2. There is no restriction on the value of the phase 
\e±{n + 1) — e±{n)\t/h ^ ngt/\/N^ provided gt/ (2A^^/^) -C 1. In fact, such phases are 
responsible for the finite corrections calculated below. 

A. Coherent cavity field 

For an initial state in which the cavity field is in a coherent state and Sz — —N/2, one 
has 

\-q,{t - 0)) = e-"'e^"«^|0) = e-^'e"(r^+^^)|0) , 
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with (5 = 0; /\/2 and 



Note that 



(5^) = ViV [(6e^^*) + (6^e-^'"*)] = 0. 
In order to compute the squeezing, we need the following averages: 



(rtrt) = Q,2gi2.tg-d2gi(2Ai-3A2)t j2 



ni 



(rt) 



^ Q,2gi2a;tgi(2Ai-3A2)tga2(e-2i^2t_i). 
^ ~2gi2a;tg-i(2Ai-3A2)tga2(e2*^2«_i). 



The value of (S) remains equal to N/2, with corrections of order a'^/N, and the squeezing, 
~ {2/\/N)^Sx is calculated as 



{ 



1 + 



COS {{2\fNg - X2)t - {ay 2) sin (2A2i)) - e" 



-2a2 sm2(A2«/2) 



+ 1 _ e-2a2sin2(A2t/2) ^^^g (|2ViV^i - sin (Aat))] j'^' 



(23) 



This expression agrees with Eq. (17) in the hmit that a <^ 1; however, it extends that result 
to all values of a for which condition (20) remains valid and for which gt/ (2A^^/^) <^ 1 
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FIG. 7: Spin squeezing as a function of gt for a = 2 and N = 60. 

[9]. Persico and Vetri [7] employ a somewhat different approach in solving this problem 
using the Holstein-Primakoff transformation and obtain a validity range, gt < ^jNjo?. 
Since 2N^I'^ > ^/Njo?, the validity range for Eq. (23) should be much greater than that 
of Persico and Vetri. To test this hypothesis, we compared the term of order in the 
exact solution with the term of (23). The two results agreed for times gtj (2A^^/^) <^ 1, 
as expected. It might be noted that Eq. (23), agrees with the exact result to order a^, 
independent of gt^ provided iV ^ 1. This is why we had to compare the oi^ terms. 

For « ^ 1, there is a slow modulation having period gt = 47r\/iV, in addition to the rapid 
oscillations having period gt — n/y/N. With increasing a, and N ^ a^, the overall period 
is gt — Att^/N, with a subharmonic having period gt — 27ry/N, and the rapid oscillations 
having period gt — n/y/N. These features are seen clearly in Fig. VA, drawn for a = 2 and 
N^60. 

Similar curves were obtained by Kozicrowski and Chumakov [5] for the field squeezing. 
With increasing a, the maximum squeezing increases slowly as is shown in Fig. VA, where 
the condition N ^ is maintained as a is varied. 

In contrast to the a <^ 1 case, the optimal squeezing for a ^ 1, always occurs at a time 
gt \fN jo?!'^ ^ \fN. In the hmit that a » 1 and z = a^gt/2VN ^ ^/a, one can show 
that Eq. (23) can be approximated as 

+ z sin {az - z) + z^sin^ [{az - z) /2]}^^^ , 

where a — AN/a'^. From this expression it is possible to show that the squeezing parameter 
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FIG. 8: Optimal spin squeezing (Cx)min ^ ^ function of aiox N ^ o? . 

goes to zero with increasing a, but that the approach to zero is slower than or^l"^ (the actual 
dependence seems to be close to Even though the field is getting more classical 

with increasing a, quantum fluctuations in the field still lead to increased squeezing with 
increasing a. Of course, if we explore the range > N, we would find a decrease in squeezing 
with increasing a, as we found for the case N — 2. 

VI. SUMMARY 

It has been shown that a linear interaction Hamiltonian between a coherent state cavity 
field and an ensemble of two-level atoms can produce spin squeezing. Analytical solutions 
for small values of the amplitude of the field state were derived, showing a reduction is 
the squeezing parameter quadratic in a. Computer simulations were used to find the best 
value for squeezing, when a is varied over a range of real, positive values. The limit of 
a large number of atoms was also examined. For an initial coherent state for the cavity 
field, it was found that the squeezing approaches zero with increasing a. This might seem 
like a remarkable result since the coherent state closely resembles a classical field for large 
a. Even though a is large, the number of atoms is assumed to be much larger than a^; 
as such the field can be totally depleted. The entanglement of the field and the spins can 
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produce significant phase shifts that can lead to spin squeezing. Although approaches 
zero with increasing a, the ratio — ^x/V^ that relates the squeezing to the Heisenberg 
limit, decreases with increasing a. If squeezing relative to the Heisenberg limit is used as a 
measure, the best squeezing is obtained for N = 2. This is in marked contrast to the optimal 
squeezing that can be obtained with nonlinear spin interactions [1, 2]. 

The interaction with a squeezed cavity field was also investigated. While a squeezed vac- 
uum field has the potential to transfer significant spin squeezing to the atoms, the degree of 
spin squeezing produced is still well above the Heisenberg limit. By constructing alternative 
squeezed states, we were able to improve the squeezing relative to that of a spin-squeezed 
vacuum, but the ultimate degree of spin squeezing that can be transferred to the atoms via 
an interaction with a cavity field remains an open question. 
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Figure Captions 

Fig. 1. Spin squeezing ^-^ as a function of gt for a — 0.4 and N — 2. 

Fig. 2. Optimal spin squeezing {^x)uiin ^ ^ function of a for = 2. The time range 
out to gt — 5000 was explored in obtaining the minimal squeezing. In this and other plots, 
the point represent actual values for which the squeezing was calculated. A line is drawn 

through these points. 

Fig. 3. Optimal field squeezing {C,q)^i^ as a function of a for = 2. The time range out 
to gt = 5000 was explored in obtaining the minimal squeezing. 

Fig. 4. Optimal spin squeezing (C2/)niin ^ ^ function of the squeezing parameter r for an 
initially squeezed cavity field for = 2. The time range out to gt — 5000 was explored in 
obtaining the minimal squeezing. 

Fig. 5. Optimal spin squeezing (^x)inin as a function of for a = 0.5. 

Fig. 6. Optimal spin squeezing (Ca:)i„in ^ function of a for A^ = 20 and < ^^t < 10. 
Since only a restricted range of gt was considered, the values plotted may not represent the 
global optimal squeezing, but still reflect the qualitative variation of {Cx)nim with a. 

Fig 7. Spin squeezing ^-^ as a function of gt for a = 2 and N — 60. 

Fig. 8. Optimal spin squeezing {^x)„iin ^ ^ function of a ior N o? . 
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